SIMILAR SUBLATTICES OF THE ROOT LATTICE 



MICHAEL BAAKE, MANUELA HEUER, AND ROBERT V. MOODY 

Abstract. Similar sublattices of the root lattice A4, are possible ^ for each index that is 
the square of a non-zero integer of the form + mn — . Here, we add a constructive 
approach, based on the arithmetic of the quaternion algebra HI(Q(-\/5)) and the existence of 
a particular involution of the second kind, which also provides the actual sublattices and the 
number of different solutions for a given index. The corresponding Dirichlet series generating 
function is closely related to the zeta function of the icosian ring. 



1. Introduction 

Any lattice F in Euclidean space possesses sublattices that are images of F under a 
non-zero similarity, i.e., under a linear map S of R*^ with {Su\Sv) = c{u\v) for all u,v ^ M*^, 
where c > and (.|.) denotes the standard Euclidean scalar product. If F' = SF C F, F' 
is a similar sublattice (or similarity sublattice, SSL) of F, with index [F : F'] = d^l"^ . If the 
lattice F has a rich point symmetry structure, there are usually interesting SSLs, beyond the 
trivial ones of the form mF with m G N. Several examples have been investigated in the 
past, compare [SI O [9l 0] and references given there. Of particular interest are root lattices, 
due to their interesting geometry and their abundance in mathematics and its applications. 
In [9], the possible values of c were classified for many root lattices, while more complete 
information in terms of generating functions is available in dimension (i < 4, see OH]. 

One important example that has not yet been solved completely is the root lattice A/^. It 
is the purpose of this short article to close this gap. After a short summary of a more general 
nature, we review the mathematical setting in 4-space via quaternions. This permits the 
solution of the problem in Section HJ using the arithmetic of a maximal order of the division 
algebra BI(Q(-v/5))- Related results will briefly be summarized afterwards, followed by some 
details on the underlying algebraic structure that we found along the way. 

2. Generalities 

A lattice F in Euclidean space M'^ is a free Z-module of rank d whose R-span is all of M*^. 
Two lattices in are called similar if they are related by a non-zero similarity. For a given 
lattice F C M*^, the similar sublattices (SSLs) form an important class, and one is interested 
in the possible indices and the arithmetic function gp{m) that counts the distinct SSLs of 
index m. Clearly, ^^(1) = 1. One can show [12j that gp{m) is super-multiplicative, i.e., 
gp{mn) > gp{m)gp{n) for coprime m,n. It is not always multiplicative, as one can see from 
considering suitable rectangular lattices in the plane [22j, such as the one spanned by (2,0) 
and (0,3). However, in many relevant cases, and in all cases that will appear below, gp{m) 
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is multiplicative. This motivates the use of a Dirichlet series generating function for it, 

grim) 



(1) Dris) := 



m=l 



which will then have an Euler product expansion as well. 

If is a generic lattice in W^, the only SSLs of F will be its integrally scaled versions, i.e., 
the sublattices of the form mr with m G N. In this case, as mr has index m*^ in F, the 
generating function simply reads 

(2) Dr{s) = (ids), 

where is Riemann's zeta function. If SSLs other than those of the form mF exist, each 
of them can again be scaled by an arbitrary natural number, so that 

(3) Dris) = Cids)D^,^{s), 

where Dp{s) is the Dirichlet series generating function for the primitive SSLs of F (we shall 
define this in more detail below). 

Let us first summarize some general properties of the generating functions in this context. 
The following statement follows from a simple conjugation argument. 

Fact 1. If F and A are similar lattices in M*^, one has Dr{s) = Da{s). □ 
Given a lattice F C M'^, its dual lattice F* is defined as 

(4) F* := {x G M"^ I {x\y) e Z for all y £ F} . 
Fact 2. If F is a lattice in M'^, one has Dr{s) = Dr*{s). 

Proof. If S is any (non-zero) linear similarity in M'^, one has 

SF C F S^F*CF*, 

which, in view of (jl]), is immediate from the relation {x\Sy) = {S^x\y). 

Observing det(S') = det(S'*), one thus obtains an index preserving bijection between the 
SSLs of F and those of F* , whence we have gp{m) = gp,{m) for all m G N. This gives 
Dr{s) = Dr*{s). □ 

Remark 1. Defining the transposed similarity via the scalar product makes it transparent 
that there is a certain scaling degree of freedom to choose the non-degenerate bilinear form 
without changing the result of Fact [21 which is in line with Fact [TJ 

3. Mathematical setting in 4-space 

As is apparent from the treatment in p|, the usual description of as a lattice in a 4- 
dimensional hyperplane of M^, see Figure [H gives access to the possible values of c, and hence 
to the possible indices, but not necessarily to the actual similarities and the SSLs. Therefore, 
we use a description in instead, which will enable us to use quaternions, based on the 
inclusion of a root system of type within one of type see jlS l Prop. 2] or ^ Thm. 4.1]. 
In particular, we shall need the Hamiltonian quaternion algebras IH(M), IHI(Q) and M.{K), 
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Figure 1. Standard basis representation of the root lattice j44. 

with the quadratic field K = Q(-v/5), generated by the basis quaternions l,i,j,k. The latter 
observe Hamilton's relations 

= f = = ijk = -1 

and will be identified with the canonical Euclidean basis in 4-space. Quaternions are thus also 
written as vectors, x = (xg, x^^, Xg), and the algebra is equipped with the usual conjugation 
X 1-^ x, where x = (xg, — x^^, — X2, — X3). For later use, we also introduce the reduced trace 
tr(x) and the reduced norm nr(x) in a quaternion algebra via 

tr(x) := X + X , nr(x) := xx = |x|^. 

If we work over the field K = Q(\/5), we shall also need its trace and norm, defined as 

Tr(a) := a + a' , N(a) := aa , 

where ' is the algebraic conjugation in K, defined by \/5 ^ — \/5- 
The first benefit of this setting is the following result. 

Fact 3. All non-zero linear similarities in ~ BI(M) can he written as either x ^ pxq 
{orientation preserving case) or as x ^ pxq {orientation reversing case), with non-zero 
quaternions p,q G EI(]R). The determinants of the linear maps defined this way are ±.\p\'^\q\'^ . 
Conversely, all mappings of the form x ^ pxq and x ^ pxq are linear similarities. 

Proof. This follows easily from the corresponding result on orthogonal matrices in M^, as every 
non-zero similarity is the product of an orthogonal transformation and a non-zero homothety, 
see [131 13] for details. The result on the determinants is standard. □ 

For some purposes, it is convenient to refer to the standard matrix representation of the 
linear map x 1— > pxq, as defined via {pxqY = M{p, g)x*. Details can be found in [T31 H]. 
Following [8], we consider the lattice 

(5) L= ((1,0,0,0), i(-l, 1,1,1), (0,-1, 0,0), i(0,l,r-l,-r)>^ 

where r = (1 -|- \/5)/2 is the golden ratio. Due to Tr(T) = 1 and N(r) = —1, r is an integer 
in K. In fact, r is a fundamental unit of Z[r], the ring of integers of K. 
The Gram matrix of L, calculated with (.|.), reads 

^2-10 0\ 
-12-10 
0-12-1 

V 0-12/ 

which shows that L is a scaled copy of the root lattice ^4 in its standard representation with 
basis vectors of squared length 2, compare [lO] and Figure [H Our choice of this particular 
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realization of the ^4 root lattice is motivated by the observation that L is a subset of the 
so-called icosian ring I, which has a powerful arithmetic structure. A concrete representation 
of I is obtained via the vectors 

(7) (±1, 0, 0, 0) , i(±l, ±1, ±1, ±1) , i(±r', ±T, 0, ±1) 

and all their even coordinate permutations. Together, they form a set of 120 vectors (all of 
length 1) which constitute a non-crystallographic root system of type H^^, denoted by Ah^^-, 
see [51 Hj for details. The orthogonal transformations {x ^ pxq , x ^ pxq \ p,q € form 
the corresponding Coxeter group Wh^ of all (120)^ point symmetries of Z\_f/4- 

The icosian ring I is the Z-span of Ah^, I := {Ah4)z, and has rank 8. At the same time, I 
is a Z[r]-module of rank 4, and can alternatively be written as 

(8) 1= ((1,0,0,0), (0,1, 0,0), i(l,l,l,l),i(l-r,r,0,l)>^[^]. 

In this formulation, I is a maximal order in the quaternion algebra M{K), which has class 
number 1, see [21]. Consequently, all ideals of I are principal, and one has a powerful notion 
of prime factorisation. 

Remark 2. The quadratic form tr(xy) = 2{x\y) takes only integer values on the icosian ring, 
and is thus the canonical form on I in the context of root systems. Relative to this form, L is 
the root lattice A4, with the weight lattice as its dual (defined relative to the trace form) 
and the index relation [^4 : A^] = 5, see [lOK^ for details. 

The icosian ring I also contains the Z[r] -module 

(9) ^ := ^ Z[Tf 

as a submodule of index 16. The same statement applies to F, which is another maximal 
order in M{K). Note that C is an index 4 submodule of I n F, the latter thus being of index 
4 both in I and in F. 

The benefit of this approach will be that we can use the arithmetic of the ring I. Let us 
recall the following result from [3], which is also immediate from Eqs. ([5]) and Q. 

Fact 4. 21 C £, i.e., x implies 2x G £. □ 

The golden key for solving the SSL problem is the mapping 7 : El(i^r) — > M.{K), x ^ x 
(also denoted by r/ later on), with 

(10) X := (xq, X]^, X3, 3:2) . 

The relevance of this map was noticed in |12j as a result of explicit calculations around 
SSLs with small indices. The rather strange looking combination of a permutation of two 
coordinates with algebraic conjugation of all coordinates, which we shall call a twist map 
from now on, is an algebra involution of the second kind [14| Ch. 1.2]. It has the following 
important properties (see also Section [6] for more). 

Lemma 1. The twist map~ of (jlOp is a Q^-linear and K-semilinear involutory algebra anti- 
automorphism, i.e., for arbitrary x,y £ M{K) and a G K, it satisfies: 
(a) X + y = X + y and ax = a'x, 
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(b) xy = y X and x = x, 

(c) X = x and thus, for x ^ 0, also (x)^^ = x^^. 

It maps K , the centre of the algebra IH(i^), onto itself, hut fixes only the elements of Q within 
K. Moreover, the icosian ring is mapped onto itself, 1 = 1, as is its algebraic conjugate, I' . 

Proof. Most of these properties are immediate from the definition, and (b) can be proved by 
checking the action of the twist map on the basis quaternions l,i,j,k. The statements on K, 
I and F follow easily from the definition in dS]). □ 

Seen as a vector space over Q, Il(K) has dimension 8 and can be split into a direct sum, 
M{K) = V++V-, where 

V± := {x £ M{K) I x = ±x} 

are the eigenspaces under the twist map, with V"+ n VL = {0} and = V^V-^. Concretely, 
one can choose the 4 basis vectors of L from ([5]) also as a Q -basis for V+. 

Remark 3. Defining u = v = j — k and w = \/5(j + k), one checks that = —5, 

v'^ = —2 and w = uv = —vu, so that D := {1,u,v,w)q is a quaternion algebra over Q, of 
type ( ~^^~^ ) in the terminology of [TH Sec. 57]. Extending this algebra by adjoining \/5 
leads to M{K), so M{K) = Q{V^)^q D. The benefit [H] is that the involution 7 looks more 
canonical in this setting. In fact, u = —u = u, and analogously for v and w. Consequently, 
for a € K and x € D, one has (ax) = a'x, which clearly shows the nature of the involution. 
Note, however, that L is a subset of V+ and is thus not contained in D. 

These observations, together with the nature of I as a maximal order in M(K), suggest 
that L might actually be the Z-module of fixed points of the twist map inside I. 

Proposition 1. One has I H V+ = {x £l\x = x} = L, with L as in 

Proof. The basis vectors of L in ([5]) are fixed under the twist map, so 5; = x is clear for all 
X G L. We also know from Lemma [1] that ][ = ][. To prove our claim, we have to show that 
no element of I \ L is fixed under the twist map. 
Let L[t] be the Z[T]-span of L. It is clear that 

L[t] = L + tL = L + tL, 

where the latter equality follows from LHtL = {0}, and also that L[t] = L[t]. Observe that 
|deti<'(][)| = 1/4 and |detK(L[T])| = \/5/4, where the iiT-determinants are just the ordinary 
determinants of the bases given in ([5]) and ([8]). Using the connection of the corresponding 
index with the ratio of the Q -determinants, see |6j for details, one concludes that 

(11) [1:L[t]] = |N(V5)| = 5, 

so that I/L[t] ~ C5. 

Consequently, the twist map, which is an involution, induces an automorphism on the cyclic 
group C5. The order of this automorphism must divide 2. Since ^ v = ^ — ]i £ I \ L[t], but 
satisfies v = —v, the induced automorphism on C5 cannot be the identity. This leaves only 
inversion (i.e., k 1— > —k mod 5), which has no fixed point in C5 other than 0, so that all fixed 
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points of the twist map inside I must lie in L[t]. It is easy to see that a quaternion of the 
form X + ry with x,y £ L is fixed if and only if y = 0, which completes the argument. □ 

Remark 4. An alternative way to prove Proposition [T] would be to use Remark [5] and to 
show that all potential x £l with x = x are elements of the dual of L (taken with respect to 
tr(xy)), but that the coset representative of A\jA^ cannot be chosen in I. In fact, extending 
Aj^ and to lattices over Z[r], compare [19j . one sees that I is intermediate between them, 
but is still integral for tr(xy). When viewed as a Z[T]-module, I has class number 1, see 
|19l Thm. 3.4]. Moreover, it is a principal ideal domain and a maximal order in E[(i^). The 
intimate relation with quaternion arithmetic is the entry point for our further analysis. 

In fact, in the setting of Remarks [2] and HJ we have L = A4 and thus L[t] = Z[r] 0^ A4. 
Then, its dual over Z[t] (relative to ti(xy) again) is Z[r]0zA|, which has K-index 5 over L[t], 
where the K-index (up to units in Z[r]) is given by the determinant of the linear mapping 
that sends a basis of L[t] to one of ^[t]*. The various inclusions can now be summarized as 

2L[r] L[r] t I t L[r]* 

together with 

2L[t] C 21 C £ C I C L[t]*, 
where the integers on top of the inclusion symbols denote the corresponding ii'-indices. The 
indices over Q are then just the squares of these numbers. 
Also, as a result of Proposition [H we get that 

(12) Aa, = v+r\AH, 

with the root system as described above. Explicitly, observing r' = 1 — r, this leaves us 
with the 20 roots of the root system Aa^, 

(13) ±(1,0,0,0), ±(0,1,0,0), ±i(±l, ±1,1,1), ±i(0,±l,r-l,-r), ±i(±l, 0, -r, r-l), 

which will become important later on, see Prop. 2] and [8j for the relation between the 
corresponding reflection groups. 

4. Similarities for A^ via quaternions 

We are interested in the SSLs of the root lattice yl4. They are in one-to-one relation to 
those of the lattice L defined in Eq. ([5]). For convenience, we define IH(ir)* = El(i^) \ {0}, 
and analogously for other rings. 

Lemma 2. All SSLs of the lattice L, defined as in Eq. (jS]), are images of L under orientation 
preserving mappings of the form x 1— > pxq, with p,q M.{K)' and K = Q(\/5 )• 

Proof. As L is invariant under the conjugation map x 1— > x, we need not consider orientation 
reversing similarities. By Fact [3l all SSLs are thus images of L under mappings x ^ pxq with 
e ]H(M)V 

A simple argument with the matrix entries of the linear transformations defined by the 
mappings x ^ pxq, compare |4j, shows that all products piqj must be in the quadratic field 
K = Q(\/5 ), which leaves the choice to take p,q £ El(i^)*. □ 



SIMILAR SUBLATTICES OF A4, 



7 



We need to select an appropriate subset of mappings that reach all SSLs of L. A first step 
is provided by the following observation. 

Lemma 3. If p £ I, pLp is an SSL of L. 

Proof. By Fact [3l pLp is similar to L, so it remains to be shown that pLp C L. Observe 
that p £ I implies p £ I, so that pLp C I is clear. If x is any point of L, we have a; = a; by 
Proposition [H Using the properties of the twist map from Lemma [U one gets 

(pxp) = pxp = pxp . 

Consequently, again by Proposition [H pxp £ L, and hence pLp C L, as claimed. □ 
Proposition 2. If pLq C L with p,q £ IHI(i^)*, there is an a £ Q such that 

q = ap. 

Proof. If p,q £ M.{K)' , the inclusion pLq C L implies pxq = pxq = qxp for all x £ L, hence 
also xpq~^ = {q)^^px = {pq~^) x. Since 1 £ L, this implies pq~^ = {pq~^) , and we get 

xpq~^ = pq~~^x, 

still for all x £ L. Noting that {L)k = M.{K), the previous equation implies that pq~^ must 
be central, i.e., an element of K. Consequently, q = ap for some a £ K. 

Since p £ M{K), we can choose some / /? G Z[t] such that w = I3p £ I. Observing 
that (Pp) = f3'p, one sees that a pxp = j^^y wxw, where 7^ N(/3) £ Z. As wLw C L by 
Lemma [3l and since LDtL = {0}, the original relation pLq C L now implies £ Q, hence 
also a G Q . □ 

An element g G I is called I-primitive when all Z[T]-divisors of q are units. Another way 
to phrase this is to say that the I-content of q, 

(14) conti(g) := lcm{Q! £ Z[t]' \ q £ al} , 

is a unit in Z[r]. Whenever the context is clear, we shall simply speak of primitivity and 
content of q. Note that the content is unique up to units of Z[r], as Z[r] is Euclidean and 
hence also a principal ideal domain. More generally, one has to formulate the content via a 
fractional ideal, see [4] for details. 

Corollary 1. All SSLs of the lattice L are images of L under mappings of the form x f-^ apxp 
with p £l primitive and a G Q . 

Proof. By Proposition [2l we know that mappings of the form x ^ f3qxq with q £ Il{K) and 
/3 G Q suffice. Since all coordinates of q are in K = Q{^/5), there is a natural number m such 
that p := mq £ C, with £ as in Q. Then, with a := one has apxp = j3qxq whence 

the mappings x 1— > apxp and x ^ f3qxq are equal. While a is still in Q, p and p are now 
elements of £ C I. 

If p is primitive in I, we are done. If not, we know that p/c is a primitive element of I when 
c = contjd?). Simultaneously, we have {p/c) = p/c' . Since c G Z[r], we know that cc' £ Z, 
so that this factor can be absorbed into a. □ 
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At this stage, we recollect an important property of the icosian ring from [3]. 

Fact 5. Let p,q £ Il{K)' be such that piq C I. // one of p or q is an I-primitive icosian, the 
other must be in I as well. 

Proof. This follows from [31 Prop. 1 and Remark 1], where this is shown for any maximal 
order of class number one. In particular, it applies to I. □ 

As a result of independent interest, we note the following. 

Fact 6. The linear map T defined by x ^ apxp has trace aN(tr(p)) and determinant 
a^N(|p|^), and its characteristic polynomial reads 

X'^ - trace(r)X^ + AX'^ - BX + det(r) 

where A = (Tr((tr(p))^(nr(p))') — 2N(nr(p))) and B = a'^ N(tr(p) nr(p)) . 

Proof. This is a straight-forward calculation, e.g., with the matrix representation from [131 14]. 
taking into account that nr(p) = (nr(p))' and expressing the coefficients in terms of traces 
and norms. □ 

In view of Corollary [H we now need to understand how an SSL of L of the form pLp with an 
I-primitive quaternion relates to the primitive sublattices of L. Recall that A\s a sublattice 
of r if and only if there is a non-singular integer matrix Z that maps a basis of to a basis 
of yl, written as Ba = BpZ in terms of basis matrices, compare [7]. The index is then 

[r-.A] = |det(Z)| 

which does not depend on the actual choice of the lattice basis. A sublattice yl of i"' C M'^ is 
called F-primitive (or simply primitive) when 

(15) contj-(yl) := lcm{m G N | yl C mr} = 1 . 

This is equivalent to saying that gcd(Z) := gcdjZjj \ 1 < i, j < d} = 1. 

Proposition 3. If p Gl is I-primitive, pLp is an L-primitive sublattice of L. 

Proof. By Lemma [3l we know that pLp C L. We thus have to show that ^pLp C L implies 
m = 1. Note first that ^pLp <Z L also implies ^pL[t]p C L[t], with L[t] = L+tL as before. 
Prom (|lip . we know that 51 C L[t], so that we get 

^P^P C ^pL[t]p C L[t] C I. 

Since p is I-primitive by assumption, then so is p. By PactO this forces 5/m to be an element 
of Z[t]. With m G N, this only leaves m = 1 or m = 5. 

Observe next that 2L[t] C C = Z[r]'*. On the other hand, it is easy to check explicitly 
(e.g., by means of the basis vectors) that \/^C C L[t]. Together with 21 C £ from Pact [H 
this gives 

4^ pip c ^pCp C ^pL[t]p C 2L[t] C £ C I. 

Invoking Pact [5] again, we see that G ^[t], which (with m G N) is only possible for m|4. 
In combination with the previous restriction, this implies m = 1. □ 
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Combining the results of Corollary [1] and Proposition [3l we obtain the following important 
observation. 

Corollary 2. The L-primitive SSLs of L are precisely the ones of the form pLp with p an 
I-primitive icosian. 

Proof. After Proposition [3l it remains to show that every primitive SSL M of L is of the form 
pLp for some primitive p G I. By Corollary [H M = apLp with a G Q and p £ I primitive. 
By Proposition [3l pLp is already a primitive sublattice, so q = ±1. □ 

Next, we need to find a suitable bijection that permits us to count the primitive SSLs of L 
of a given index. Recall from |15^ [8] that the unit group of I has the form 

r = {x e I I N(nr(3;)) = ±1} = {±r™e | m G Z and e E AhJ , 

where Ah^ is defined as above via Eq. ([7]) and satisfies Ah^ = n S'^ = I n S'^, see [U [15] for 
detail^. 

Lemma 4. For p £l, one has pLp = L if and only if p G I^. 

Proof We have pLp C L for all p G I by Lemma El The corresponding index is given by the 
determinant of the mapping x 1-^ pxp. By Fact [6l the equality L = pLp is then equivalent to 
N(IpI^) = 1 ^-iid hence to N(nr(j))) = ±1, which, in turn, is equivalent to p G I^. □ 

We need one further result to construct a bijective correspondence between primitive SSLs 
of L and certain right ideals of I, which will then solve our problem. 

Lemma 5. For primitive r, s G I, one has rl = si if and only if rLr = sLs. 

Proof. Since r, s G I, it is clear that rl = si I = r^^sl =^ r~^s G I. Similarly, s~^r G I, 
so r^^s G I^. Lemma m now implies r~^sL{r~^s) = L, which gives rLr = sLs. 

In the reverse direction, suppose that rLr = sL's, which gives yLy = L with y := r~^s. 
Choose a & K so that ay is a primitive element of I, which is always possible. Then, 

ayLay = aa'yLy = aa'L 

is a primitive sublattice of L by Proposition [3l whence aa' = ±1 and ayLay = L. 
This implies e := ay = ar~^s G I^ by Lemma [H Then, 

r = ree~^ = ase~^ G al, 

so that a'r G I due to aa' = ±1, where a' £ K hj construction. Since r G I is primitive as 
well, such a relation is only possible with a G Z[t], in view of the properties of the I-content 
of r. Consequently, aa' = ±1 now gives a G Z[t]^, so that y is an element of I. Lemma [J] 
now implies y G , whence yl = I, and finally si = rl. □ 

In view of our discussion so far, we call a right ideal of I primitive if it is of the form pi for 
some primitive p £l. 

^Let us take this opportunity to mention that [4], in deviation from other conventions and from the one 
adopted here, uses the notation I*^ for the I-units on S^, i.e., for the elements of Ah4 only. 



10 



MICHAEL BAAKE, MANUELA HEUER, AND ROBERT V. MOODY 



Proposition 4. There is a bijective correspondence between the primitive right ideals of I 
and the primitive SSLs of L, defined by pi ^ pLp. Furthermore, one has the index formula 

[I : pi] = N(nr(p)2) = lSS{\p\^) = [L : pLp] 

under this correspondence. 

Proof. It is clear from Lemma O that the mapping is well-defined and injective, while Corol- 
lary [2] implies its surjectivity. 

The index of pi in I derives from the determinant formula in Fact [3] by taking the norm in 
Z[r]. The index relation now follows from Fact[6j □ 

Remark 5. The results of Lemmas U] and O also show that zAjy^ = R which is a subgroup 
of I^ of order 120, can be viewed as the standard double cover of the rotation symmetry group 
of L, which has order 60 and contains the point refiection in the origin (note that —1 is a 
rotation in 4-space). This is a geometric interpretation of the bijective correspondence. Let 
us, in this context, also point out that pIpD V+ = pLp, but that pi n V+ is generally a much 
bigger set. 

With the representation of a general SSL as an integer multiple of a primitive SSL men- 
tioned earlier, and observing that all possible indices of SSLs are squares, we can now solve 
our original problem. We need the Dirichlet character 



(16) x{n) 



0, n = (5) , 

1, n = ±l{5), 
-1, n = ±2(5). 

Note that the corresponding L-series, L(s,x) = Ylm=i xii^)^~'^j defines an entire function 
on the complex plane. Furthermore, we need the zeta functions of the icosian ring, 

(17) (lis) = Ci^(2s)Cx(25-l) and C,j{s) = Ci^(4s) , 

where Ci(s) and Cii('5) denote the Dirchlet series for the one-sided and the two-sided ideals 
of I, respectively. Here, Cxi^) ~ C{s) L{s,x) is the Dedekind zeta function of the quadratic 
field iC, see [211 H E] for details. 

Theorem 1. The number of SSLs of a given index is the same for the lattices A4 and L. 
The possible indices are the squares of non-zero integers of the form k'^ + k£ — £'^ = N{k + iT). 
Moreover, all possible indices are realized. 

Furthermore, there is an index preserving bijection between the primitive SSLs of A4, and 
the primitive right ideals of I. When f{m) denotes the number of SSLs of index m? and 
/P''(m) the number of primitive ones, the corresponding Dirichlet series generating functions 
read 

J->aAs) ■= 2^ — oT ~ C(4gj - — ; r 



m=l 



and 



m=l 
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with Ck('^) '^"'^ -^(^' '^^ defined above. In particular, each possible index is thus also realized 
by a primitive SSL. 

Proof. The first claim follows from Fact [H while the index characterization follows either 
from [9] or from the index formula in Proposition [H Recall that k'^ + ki — i'^ = N(A; + ir) 
is the norm the principal ideal {k + £r)Z[r], and that all Z[r]-ideals are of this form. Since 
the reduced norm as a mapping from icosian right ideals to ideals in Z[t] is surjective, each 
possible index is realized. 

In view of the previous remarks, we have 

where, due to our bijective correspondence and the index formula of PropositionHJ the various 
formulas for the Dirichlet series are immediate, see [H |6] for a derivation of C]f'^('5) as the 
quotient stated above. The last equality then follows from (|17p . Since f^'^{m) vanishes 
precisely when f{m) does (see below for an explicit formula), the last claim is clear. □ 

Inserting the Euler products of C{s) and Cit(s)) one finds the expansion of the Dirichlet 
series Da4,{s) as an Euler product, 

1 -TT l+p-'^' 1 T-r l+p-'^' 1 



(1-5~2^)(1 



n 



^ p=±l (5) 



P 



-2s 



{l-p 



l-2s^2 



n 



1 - p-^' 1 - p 



2-4s 



' p=±2(5) 

and similarly for Z)^'^(s). Consequently, the arithmetic function f{m) (and also /f"'(m)) is 
multiplicative, i.e., f{mn) = f{m)f{n) for m,n coprime, with /(I) = 1. The function / is 
then completely specified by its values at prime powers p^ with r > 1. These are given by 



/(/ 



2 (l-p''+i)-(r+l)(l-p^)p'' 
2-p''-p''+^ 



1-p^ 



if p = 5, 
for primes p = ±1 (5), 
for primes p = ±2 (5) and r even, 
for primes p = ±2 (5) and r odd. 



and similarly for f'^^(m): 

'6- 5'^-^ 



rip' 



(r + l)p^ + 2rp^ + [r — l)p 
0, 



r-2 



if p = 5, 

forp = ±1 (5), 

for p = ±2 (5) and r even, 

for p = ±2 (5) and r odd. 



The first few terms of the Dirichlet series thus read 



1 + 
1 + 



+ 
+ 



I 11 I 24 I 26 I _40 I 36 I 31 I 
-I- g2i- -r -I- -I- Yg^J -I- "T 25^ 29' 



60 



66 



_ I 10 I 24 I 

where the denominators are the squares of the integers previously identified in [9], see also 
entry A031363 of [20]. Further details can be found in |12] . 

These series may now be compared with the zeta functions Ci(s) and Cj^s), which reveals 
the origin of the various contributions. In particular, the 6 SSLs of index 16 stem from the 



20 



^ Tg^ ^ 20 



30 



+ 



30 



+ 



60 



+ 



64 



+ 



50 



+ ... 
+ ... 
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5 generators of primitive ideals of I of index 16 together with the SSL 2^4. No such extra 
solution exists for index 25, while index 81 emerges also from the SSL 3^4. 

Due to our definition with f{m) being the number of SSLs of A4 of index m^, the Dirichlet 
series generating function of the arithmetic function / is Da^{s/2), which has nice analytic 
properties. In particular, it is analytic on the half-plane {a > 2}, where we write s = a + it as 
usual. Moreover, it is analytic on the line {fi = 2}, except at s = 2, where we have a simple 
pole as the right-most singularity of DAi{s/2). Consequently, one can derive the asymptotic 
growth of /(m) from it. Since the value of the arithmetic function /(m) fluctuates heavily, 
this is done via the corresponding summatory function 

:= ^ f{m) ~ e'y, as x ^ 00, 

m,<x 

where the growth constant is given by 

g = Tess=2 = ^^^^^^^^ = ^Vb logir) ~ 0.538011. 

This result relies on Delange's theorem, see [U Appendix] for details. The corresponding 
calculations for the asymptotic behaviour of /P'^(m) are analogous. 

5. Results for related lattices 

From previously published results, one can read off or easily derive the generating functions 
for the root lattices with d < 3. 

Theorem 2. The Dirichlet series generating functions for the number of SSLs of the root 
lattices with d < 3 are D^_^{s) = C(s), Dj^^{s) = Cq(^3); ''^ith ,^3 = e^'^*/^ and Cq(^^^-^{s) 
Dedekind zeta function of the cyclotomic field Q(^3), and Dy^^{s) = C(3s) <Pcub(3s), where 

is the generating function of the related cubic coincidence site lattice problem derived in [UlG]. 

Proof. Ai is a scaled version of the integer lattice Z, whence Dai{s) = Dz{s) = C{s) follows 
from Fact [1] together with [3, Prop. 3.1]. The triangular lattice A2 is a scaled version of the 
ring of Eisenstein integers in Q(i^3), so that this generating function follows from [21 Prop. 1]. 
It can also be written as a product of two Dirichlet series, 

one for the scaling by integers and the other for the primitive SSLs. 

For the primitive cubic lattice Z^, the SSL generating function was previously identified as 
C(3s) ^cub(3s) in [21 Thm. 5.1]. Noting that A3 is a version of the face-centred cubic lattice, 
compare [10], one has to check that the cubic lattices share the same SSL statistics. This is a 
straight-forward calculation with their basis matrices and the integrality conditions for simi- 
larity transforms, similar to the one outlined in |3i Sec. 5], using the Cayley parametrization 
of matrices in S0(3,Q) from [T], see also [12] for details. □ 
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Clearly, by Fact [2l one then also has the relation 
(18) Da^{s) = Z)a-(s) = D^^is). 

This identity can also be derived [22] by a consideration of the shelling structure of A'^ relative 
to Z^, which bypasses the explicit basis calculations referred to above. 

Other completely worked out examples include the square lattice Z^, with Dj2{s) = 
Cq)(^j)(s), see [3], and the cubic lattices in d = 4, see [1]. Also, various Z- modules of rank 
r > d in dimensions d < 4 are solved in O HI O [2] . Common to all these examples is the 
rather explicit use of methods from algebraic number theory or quaternions in conjunction 
with suitable parametrizations of the rotations involved. In higher dimensions, results are 
sparse, compare [9j and references given there, and we are not aware of any complete solution 
in terms of generating functions at present. 

6. Algebraic afterthoughts 

The special twist map T played a crucial role in solving the sublattice problem above. It 
is an odd looking map at first sight, but has a number of interesting properties that are 
worth considering in more detail. To do so, we also write r}{x) = x from now on, and use 
both notations in parallel. Moreover, we define inner automorphisms Ta of the quaternion 
algebra El(i^) via Ta{x) = axa~^, for a £ E[(i^)*. Finally, we call any iiT-semilinear involutory 
anti- automorphism of M{K) a twist map. 

Lemma 6. If (T^ and are two twist maps of M(K) that map I into itself, their product is 
an inner automorphism of I, i.e., o-^a2 = for some e G I^. In particular, any such twist 
map is of the form T^ij for some e G I^. 

Proof. Since both ai are if-semilinear, their product is i^-linear. Also, the product cT^Cg of 
two anti-automorphisms is an automorphism of IHI(i^), which must be inner by the Skolem- 
Noether theorem [141 Thm. 1.1.4]. This gives o"^(T2 = Ta with a £ M{K)'. By construction, 
Ta is also an automorphism of I, whence we may assume a G I without loss of generality. 
This implies al = la, so that al is a two-sided ideal. Within I, this implies a = ae with 

7^ a G Z[t] and e G I^. Observing T^e = establishes the claim. □ 

For our further arguments, it is advantageous to consider the mapping E[(Er) — > V+, 
defined by x 1-^ xr]{x) = xx, where 1+ = {x G M.{K) | x = x} as before. Clearly, one has 
C L, as a consequence of Proposition [H 

Lemma 7. The mapping satisfies = Tl){Ah^) C Aj^^. Moreover, the corresponding 

restriction, -d/^: Ah^ — > ^A^, is a surjective 6-to-l mapping. 

Proof. Any a G has the form a = ztr^e for some m G Z and e G Ah^, so that aa = 
[tt')'^ E£ = (— l)™ee. The identity = '^{Ah^) now follows from the observation that 

1 G Ah4^, with i?(i) = i^ = —1. If e G Ah^, ££ lies both in L and in Au^, and thus in 
Aa^ =LnAH^, see Eq. 

Observe next that the equation ee = 1 has precisely 6 solutions with e G Ah^, namely 
±1 together with ±^(±1, 0, r, r — 1). Under multiplication, they form the cyclic group Cg, 
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generated by the quaternion ^(l,0,r,r — 1). Clearly, one has ^{a) = 'd{b) with a,b £ Ah_^ if 
and only if 6 = ae with e G Cg, so that a coset decomposition of the 120 roots of with 
respect to Cq establishes the second claim. □ 

Define now, for an arbitrary e G I^, the mapping rj^ = T^rjT^-i. Since ij^e = rj^ for any 
7 € Z[t]^, all such maps can be obtained from an e G Ah^- 

Proposition 5. The twist maps of M{K) that map I into itself are precisely the maps % = 
TsrrjT^-i = T^-Tj, with e G Ah4^. Equivalently, they are precisely the maps Taf] with a G Aa^- 

Proof. By Lemma El we know that any twist map is of the form Taf] for some a G I^. It is easy 
to check the commutation relation r]T^-i = ^a^) which also implies, via left multiplication by 
Ta, that 

A twist map is an involution (in particular, one has rj^ = Id), so that (Tarj)'^ = Ta(a)-'^ has 
to be the identity. Since T;, = Id = Ti with 6 G implies b £ Ci K = Z[r]^, the previous 
condition means a = f3a for some (3 G Z[r]^. Consequently, a = {f3a) = j3'a = P'(3a, whence 
/3' f3 = 1. This implies /3 = itr^'" for some m G Z. Observing that Tt^^ = Ta, we may 
thus, without loss of generality, replace a by T™a, so that the involution condition reduces to 
a = ±a. 

Next, write a = with 7 G Z[r]^ and e G Ah^, which is possible due to the structure of 
I^. Inserting this into a = ±a results in 7^"^7' = =bee~^. Observe next that Z[t]^ and Ah^^ 
share only ±1. Within Z[t]^, 7~^7' = —1 has no solution, while e = —s has no solution in 
Ahj^, meaning that we must have 7~^7' = 1 = ee~^. This leaves us with the condition a = a, 
which is precisely solved by the 20 elements of Aa4^, see Eq. (jl2p . All possible twist maps are 
thus of the form Tai] with a G Aa^- 

Observing ee = t?(e), the first claim now follows from Lemma [3 □ 

Theorem 3. There are exactly 10 distinct twist maps of M{K) that map I into itself, which 
may be written as % with e running through the 10 positive roots of Aa^ ■ 

Moreover, the fixed points of ije are precisely the points of the lattice T^{L), each of which 
is a root lattice of type ^4. 

Proof. The twist maps TaTj with a G Aa^ from Proposition [5] form precisely 10 distinct maps 
because Ti, = Ta holds here if and only if 6 = ita. Consequently, the restriction to the positive 
roots of Aa^ suffices. 

Prom Proposition [H we know that L = {x G I | x = x} is the lattice of fixed points of 
rj^ = rj. It is an immediate consequence that r]^ fixes the points of eLe~^ = T^i^L). □ 

Remark 6. There are 60 congruent copies of L within the icosian ring I, each of which must be 
of the form aLb with a,b £ Ah^ due to the structure of . Observing aLa = L = —L = —aLa 
for a G Ah_^, it is evident that the stabilizer of L, when seen as a subgroup of Ah^^ x Ah^, 
has order 240. Since the order of Ah^ x Ah_^ is 120^, the length of the orbit of L under the 
action of this group is 60, which are the congruent copies. Among them, the 10 lattices of 
Theorem [3] are precisely the ones that contain 1. They are thus special in the sense that they 
line up with the arithmetic of I. 
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Consider T : — > Aut(I), e ^ T^, which is a group homomorphism. As = T_£, 
T{A}{4) — y is the standard icosahedral group of order 60. If c denotes the conjugation map, 
X 1-^ c{x) := X, one has cT^ = T^c for all e € Ah^^ so that 

(r(z\H4),c> ~ y XC2 = y/. 

is the full icosahedral group of order 120. By Lemma [H c and 77 commute as well, while 
'qTf.r] = T| shows that 

{T{AH,),c,ri) ~ n X (7?) = (y X (r/)) x (c) . 

Lemma 8. One /las y x [rj) — S5, with S5 being the full permutation group of 5 elements. 

Proof. Note that Y is isomorphic with the alternating group of 5 elements, which is simple. 
As (rj) — C2, y X (r/) is a C2-extension of Y, where the automorphism on Y induced by ry is 
Te I— > Ti, hence not the identity. As the automorphism class group of Y is the cyclic group 
C2, there are (up to isomorphism) only two C2-extensions of Y, namely Y^ and 5*5. Since 
y^j ~ y X C2, the claim follows. □ 

Let us add some comments on the geometric meaning of the 10 twist maps, in the setting of 
root systems, where we focus on the interplay between twist maps that map I into itself and 
elements of order 3. Here, we use the natural quadratic form tr{xy) mentioned in Remark [2j 

The intersection of the fixed points of rj and the pure icosians Iq = {2; S I | x + x = 0} is 
the Z-span of a root system of type A2. The latter consists of the 6 roots 

±(0,1,0,0), ±i(0,±l,r-l,-T), 

which is also the intersection of Iq with our original A4 root system from Eq. (jl3p . 

Now, the twist maps are semilinear involutions of I, each of them a conjugate of our given 
one, r], by an element of the group I^. Since all conjugations by elements of are orthogonal 
transformations that stabilize 1 and Iq j each twist map likewise fixes an A2 root system which 
lies both inside its own A4 and inside Iq- 

We can easily see the geometric meaning of these 10 root systems of type A2. We have 
noted that each a € Z\//4 determines, by conjugation, an orthogonal transformation Tq that 
preserves I, Iq, and Ah^, where Ah^ = InS^ as mentioned above. In particular, Tq preserves 
lo n Ah^, which is a root system of type H3, with 30 elements, that we denote by AH3. Its 
convex hull is the icosidodecahedron, a semi-regular polytope with the 30 roots as its vertices, 
20 triangular faces, and 12 pentagonal faces. It is the dual of Kepler's famous triacontahedron, 
see im Plate I, Figs. 10 and 12] for details. The 20 triangular faces, which come in centrally 
symmetric pairs, give rise to 10 axes of 3-fold symmetry — which, in fact, correspond to the 
elements of order 3 in the symmetry group Yh of the icosahedron. We shall show that each of 
these axes has a hexagon of roots in the plane orthogonal to it. These hexagons are precisely 
the root systems of type A2 that arise from our twist maps. 

To see this, recall that the elements of the form Tq, with a G Ah4^, give rise to the symmetry 
group (of order 120) of the root system Ah^, via the mapings x Ta{x) and x 1— > Ta{x). Here, 
the first type provides the 60 orientation preserving elements (isomorphic to the alternating 
group on 5 elements), and the second type the 60 orientation reversing elements. Together, 
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they form Y^, the symmetry group of the icosahedron (and hence also of the icosidodecahedron 
and the triacontahedron) . 

Because all our twist maps are conjugate to one another, it suffices to look at the A2 lattice 
indicated above. It has the A2 root basis 

{(0, -1,0,0), i(0,l,r-l,-r)}, 

and the product of the two reflections generated by these quaternions is the rotation of order 
3 that is a Coxeter element of this root system (the only other one being its inverse, which is 
obtained by multiplying the two elements in the opposite order). It is the inner automorphism 
determined by 

-i(0,l,0,0)(0,l,T-l,-T) = i(-l,0,r,-r') =:2. 

Of course, z is an element of Ah^, and it is easy to see that = — 1, so has order 3. These 
inner automorphisms and their inverses give 20 elements of order 3, which then exhaust the 
elements of order 3 in Yh. This establishes the connection between the A2 root systems in 
Aff^ and the rotations of order 3 in the icosahedral group. 

Each A2 inside Ah-^ completes uniquely to an ^4 which contains ±1. This is easy to see 
from Eq. ([5]), in which the four icosians displayed are an ^4 basis which includes the A2 that 
we have been using in our above argument. 

So, in short, geometrically we can locate the A4 root systems that are the fixed points 
of the special twist maps as follows. Each of the great circles through 6 vertices of the 
icosidodecahedron of roots in Iq is a root system of type A2, and it extends uniquely to a root 
system of type A4 inside Ah^ that contains 1. 
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